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The Extended Analysis On New Generalized Chaplygin Gas 

WANG Jun*, WU Ya-Bo^ WANG Di, YANG Wci-Qiang 
Department of Physics, Liaoning Normal University, Dalian 116029 

We will extend the study of the new generalized Chaplygin gas (NGCG) based on [JCAP 
0601(2006)003]. Concretely, we will not only discuss the change rates of the energy densities and the 
energy transfer of this model, but also perform the Om diagnostic to differentiate the ACDM model 
from the NGCG and the GCG models. Furthermore, in order to consider the influence of dark 
^ , energy on the structure formation, we also present the evolution of growth index in this scenario 

\ with interaction. 
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Recent observations of type la supernovae (SNc Ia)[l, 2] indicate that the expansion of the Universe is accelerating 
at the present time. From the observations of large scale structure (LSS)[3], we know that the Universe is spatially 
flat. The Wilkinson Microwave Anisotropy Probe (WMAP) indicates that our Universe components are as follows: 
2 ' usual baryon matter occupies about 4%, dark matter occupies about 23% and dark energy occupies about 73%. It is 
' clear that the Universe is dominated by an exotic component with large negative pressure, referred to as dark energy. 

■ 

I— ' ' So many scientists believe the accelerating expansion of the Universe is due to the dark energy. 

fSj ' Based on this opinion many dark energy models have been proposed. The simplest form of dark energy is 

>'• . 

, cosmological constant A which would encounter "fine-tuning" problem and "coincidence" problem. Other valid dark 

I energy models are provided by scalar fields, such as: Quintessence [4, 5], which is introduced to solve the "coincidence" 

■ problem and characterized by the equation of state (EOS) Wde between -1 and -1/3 (namely, —1 < Wde < —1/3); 

• Phantom (ghost) field[6], which owns a negative kinetic energy and characterized by the EOS w^e less than -1 

C\| ' (namely, Wde < —1); Tachyon field[7, 8]which can act as a source of dark energy depending upon the form of the 

tachyon potential, and so on. Other scenarios on dark energy include brane world[9], generalized Chaplygin gas[10], 

, , , holographic dark energy [11], ect. 

^ , In the generalized Chaplygin gas (GCG) approach, dark energy and dark matter can be unified by using an exotic 

< 

- ■ ■ equation of state. This point can be easily seen from the fact that the GCG model behaves as a dust-like matter 
at early times and as a cosmological constant at late stage. Since the EOS of dark energy still can't be determined 
exactly, the observational data show the EOS of dark energy is in the range of (—1.38, — 0.86)[12, 13], the GCG model 
is generalized to accommodate any possible X-type dark energy with constant w^e- The new generalized Chaplygin 
gas (NGCG) as a dark energy model is proposed for unification of X-type dark energy and dark matter [14]. The EOS 
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of this model is p = —A{a)/p°', where a is the scale factor and A{a) = — u)rfej4a~^'^"''™''=-"^+"-'. We know that the 
NGCG model behaves as a dust-like matter at early times and as a cosmological constant at late stage. Also this 
model is a kind of the interacting XCDM model[14]. Based on the Ref.[14], in this paper we will extend the analysis 
on the NGCG model by discussing the change rates of the energy densities for dark energy (DE) and dark matter 
(DM) as well as the energy transfer. Moreover, we perform the Om diagnostic to differentiate the ACDM model from 
the NGCG and the GCG models. Furthermore, in order to consider the influence of DE on the structure formation 
in this scenario, we present the evolution of growth index in this model with interaction. 
In the framework of FRW cosmology, the EOS of the NGCG model is described as [14] 

Ma) f,. 

PNGCG — a ' 

Pngcg 

where a is a real number and A(a) is a function depending on the scale factor of the Universe. The form of the 
function A{a) is: 

i(a) = -u;de^a-3(i+--)(i+"), (2) 

where A is a positive constant and Wde, which should be taken as any possible value in the range (-1.38, -0.86), is a 
constant. We can see explicitly that when Wde = — 1, the EOS can reduce to the GCG scenario, while when a = 
the XCDM model can be obtained again. 

In the following discussions we suppose that in the NGCG model the fluid is made up of two components, one is 
the DE component marked as pde, and the other is the DM component marked as pdm, i-e. 

Pngcg ~ Pde + Pdm, (3) 

PNGCG =Pde- (4) 

Then, the energy density can be expressed as [14] 

PNGCG = PNGCGOa~^[l - ^ + ^a'^-- TT^ , (5) 

where As = -f^. The energy densities of the DE and the DM components can be respectively expressed as 

Pde = Pdeoa"'''+"'^=''+"^l [1 - ^ + ^a-3»-(i+")]T^^-i, (6) 

Pd,n = Pd,noa-''[l ~As+ ^a-3-^^(i+")]TT^-i. (7) 

Below we study the change rates of the energy densities of DE and DM. First, we consider the case without 
interaction between them. By means of the continuity equation, the energy densities of DM and DE can be respectively 
obtained as follows: 

Pdm = Pdmoi^ + z)^ , (8) 

Pde = {pNGCGoil -As+ A(l + ^)3"'-(l+")] T^J - pd„o}(l + z)^ (9) 
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Thus, the change rates of the energy densities of DE and DM can be obtained: 

dpde 



dz 



dpdm „ /I I \2 

— 1 = OpdrnO[l + Z) , 



(10) 



(11) 



where M = (1 + z)^[l - A,, + ^,,(1 + z)3"'^<=(i+")] tt^ (1 + ^^^,^^^1+!^^^°'^) ) ^nd A, = ^ = 0.74 (the current 
parameters used in this paper are: fldmo ~ 0.25, fldeo = 0.7, and ilbo ~ 0.05.). It is easy to see that the sign of the 
value of ^^4^ hes on M. Hence the evolutional trend of DE density can be discussed by the factor M as well. 




red shift z 



FIG. 1: The evolutional trends of the factor M and 



with z. Here we choose As = 0.74, a — 0.5, Wde 



1.2. 



The evolutional trends of the factor M and are illustrated in Fig.l, from which we can see that when z > 0, 
the change rate of the energy density of DE decreases faster than that of DM as the redshift z becomes low. When 
z < 0, it keeps decreasing but ^^^^ would approach to 0. 

In what follows we will consider the NGCG model with interaction between DE with constant Wde and DM. 
According to the Eqs.(6) and (7), we can give the expressions of the change rates of the energy densities of DE and 
DM as follows: 

dpde 



dz 



dp 



'dm 



dz 



SpdeoW, 



SpdmoS, 



(12) 



(13) 



where W ^ [1 + z)3-d.(i+«)+2[i _ A, + ^(1 + z)3w^M+c^)]Th~^[l + Wde{l + a) - — ^^+^^(i+^)3»,,(i+„) j, 
S = {1 + z)2[l -A,+ A,{1 + z)3-d=(i+a)]TT^-i(i _ _^^i£^(i±£i!;^l!^). From Eq.(12), it is easy to see that the 



sign of the value of depends on the factor W. Fig. 2 shows the relation between factor W and the redshift 
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FIG. 2: the relation between factor W and the redshift z. Here we choose = 0.74, a = 0.5, w^^ — —1.2. 

From Fig. 2, it is easy to see that as the redshift z becomes low the change rate of the energy density of DE is from 
increasing gradually to decreasing sharply. Specially, when z > 4 it would approach to 0. It arrives at the maximum 
when 2 = 0.53. While when z ~ 0.2 it vanishes. It follows that the point of 2 = 0.53 is a transformation point. Note 
that there is not any transformation point in the case without interaction. 

From the above discussions, we see explicitly that there would exist energy transfer between DE and DM. Now, we 
consider a feasible interaction with an ansatz F = 3Hc^p [15], where p is the total energy density and is a constant 
denoting the transfer strength. We suppose the components of the NGCG interact through the interaction term F. 
The components of the NGCG respectively satisfy 

Pde + 3H{1 + Wde)pde = "F, (14) 
Pdm + SHpdm = r. (15) 

Note that F > 0, which implies there is an energy transfer from DE to DM. By using Eqs.(7), (13) and (15) as well 
as ^ = — (1 + z)H, we can obtain: 

f = -~9HpdrnoAsWdeaB, (16) 

where B = {I + ;^)3-.e(i+o)+2[i _ + ^,(1 -f ,)3^..a+^)]Th-^[l + Wd^il + a) - ^^t!;!^!?/-^.?^" (i tT' 1 • From 
Eq.(16), we can see that the evolution of F can be described by the factor B. The relation between factor B and the 
redshift z is shown in Fig. 3. 

According to Fig. 3, we can see that when z > 0.38 the transfer direction of the energy flow is from DM to DE and 
the quantity of energy transfer is from increasing gradually to decreasing sharply as the redshift z becomes low. But 
when — 1 < 2 < 0.38 the transfer direction of the energy flow just reverses and the quantity of energy transfer is from 
increasing to decreasing sharply, until approaches to at z = —1. Clearly, the point of z = 0.38 is the transformation 
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FIG. 3: The relation between factor B and the redshift z. Here we choose = 0.74, a = 0.5, ~ —1.2. 

point and the quantity of energy transfer doesn't change any more at this point. Specially, when z = 0.78, the quantity 
of energy transfer from DM to DE reaches to the maximum. While when z = —0.02, it reaches to the maximum from 
DE to DM. It follows that the above discussions not only show the evolutionary laws of the energy transfer in the 
past and at present, but also predict the situation in the future. 

Below, we will determine the present value of c^. The ratio of the energy density can be defined as f = Thus 
by means of Eqs.(6),(7),(14) and (15) the evolution equation of ratio r can be expressed as 

( ) = 3i/r[c(- — + -^) + Wdel (17) 

Pde Pdm Pde ^ ^dni ^ ^de 

from which, the following expression is easily obtained: 



c 



,2 



-Wdeil + a) 11 



«^rfe](7^ + 7^)"^ (18) 



(1 + z)H ' rirfm Qde 

Hence wc can obtain the corresponding present values of to the different values of Wde and a shown in Tab. I. 



TOde —0.8 


-1 


-1.2 a 





0.2 


0.5 


-0.0193076 


-0.0242937 


-0.0292797 


-0.0242937 


-0.0242937 


-0.0242937 



(a)a = 0.5 {h)wde = -1 

TABLE I: The corresponding present values of c^. The current parameters are taken to be z = 0, fidmo = 0.25, fldea = 0.7, 
Ho = 70.5[13]. 



From Tab. 1(b), it is easy to see that the values of don't change with a just as the case in the ACDM model. 
Moreover, we will use a new geometrical diagnostic method, Om diagnostic [20], to differentiate the ACDM model 
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from the NGCG and the GCG models. The Om diagnostic is defined as follows 

Om{x) = — W ,x^l + z, (19) 

X'^ — \ 

where E'^{x) ^ H{x)/Ho = (1 - nbQ)x^[l - As + AsX^'"-'-^^+°'^^ + n^x^ for the NGCG model. 

For the ACDM model, Om{z) = ilom is a constant. It provides a null test of cosmology constant. The benefit of 
Om diagnostic is that the quantity Om{z) can distinguish DE models with less dependence on matter density fi„jo 
relative to the EOS of DE[20]. We can get the expression of 0'm{x) in the NGCG model as 

^ , , {\ - ni,o)x^[l - As + ^a;3^<^e(i+a)] + ^^^^^ _ 1 

Om{x) = -3—^^ (20) 

where As = 0.74. It follows that when taking a = 0.5, = —1, Eq.(20) can reduce into the Om{x) of the GCG 
model, but when taking a = O^Wde = —1, it corresponds to the Om{x) of the ACDM model. For the NGCG model, 
the curve of Om{x) corresponding to a = 0.5, Wde = —1.2 is plotted in Fig. 4. 
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FIG. 4: The Om(x) evolution diagram. The current parameters are taken to be Q.dmo ~ 0.25, fldeo = 0.7, fibo = 0.05 and 
As = 0.74. 

The Fig. 4 shows that when x > A the evolutionary trajectory of the NGCG model is similar to one of the GCG 
model, but in the recent period the both become very different. It follows that the both NGCG and GCG models are 
completely different from the ACDM model. 

Furthermore, in order to study the influence of the NGCG model on the structure formation, now we discuss the 
growth index of the NGCG model. The growth index is an important quantity to test a model, which can be measured 
by the galaxy correlation function or the peculiar velocities. Its definition is as follows: 

^ _ din (5 _ adS 
d\na S da^ 
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where S = is the matter density contrast and a is the scale factor of the Universe. 

According to the perturbation equation[16]; the growth index / satisfies the following equation: 

/'+- + (--l)/-|-^^rn=0, (22) 

a a la 

where the prime denotes the derivation with respect to the scale factor a. In the NGCG model, we take 17„i ~ ^Idm+^b- 
the corresponding expression of ^dm and fib reads [14]: 

^dm = ndrnoE-'^a-^[l -As+ ^a-3-'-(i+°)]TT^-i, (23) 

fib = ^!bo^"'a"^ (24) 
where E = {{I - nbo)a-'^[l ~ + ^a-3"'<ie(i+a)]Tir + r2toa-3}i/2 and A, = 0.74. 




FIG. 5: The evolutionary trajectories of the growth indices / for the NGCG, the GCG and the ACDM models. The current 
parameters are taken to be Q.dmo = 0.25, fideo = 0.7, Slto = 0.05 and As = 0.74. 

The evolutionary trajectory of the growth index / for the NGCG model is plotted in Fig. 5. By calculations, 
the value of / for the NGCG model is 0.592 at a = 0.87 (namely, z = 0.15). From the 2dF galaxy redshift survey 
(2dFGRS) [17-19], we know that the observational value of / is / = 0.51 ±0.1 or / = 0.58±0.11 at the effective redshift 
of the survey z ~ 0.15. It means the theoretical value of growth index given by us is consistent with observation. In 
Fig. 5, we also give the evolutionary trajectories of growth indices of the GCG model (namely, a — 0.5, Wde = — 1) and 
the ACDM model (namely, a = QTWde = ^1)- By comparing these three evolutionary trajectories in Fig. 5, we can 
find that at early times, the growth indices of these three models are the same; but near a = 0.09, the evolutionary 
trajectory of ACDM model would be different from the other two; while near a = 0.2, the evolutionary trajectory of 
the GCG model would deviate from one of the NGCG model. 
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In summary, we have extended the analysis on the new generahzed Chaplygin gas (NGCG) model as the unification 
of DM and DE. Concretely, we have studied the evolutional trends of the energy densities, Om diagnostic and the 
growth index in the NGCG model. By complicated calculations and analysis, we can give some comments as follows: 

1. In the NGCG model without interaction between DM and DE, the change rate of the energy density of DE 
decreases gradually as the redshift z becomes low, but it would be from increasing gradually to decreasing sharply in 
the case with interaction, while z = 0.53 is a transformation point. 

2. For the NGCG model with interaction term F = SHc^p, the evolutionary laws of the energy transfer have been 
discussed. Furthermore, we have also given the present values of when taking a fixed constant a (or Wde), but 
different values of Wde (or a). The results show that for a fixed Wde, the present values of don't change with a, 
which is as the same as one of the ACDM model. 

3. We have performed the Om diagnostic to the NGCG model and shown the discrimination among the NGCG, 
the GCG and the ACDM models. 

4. The evolutionary trajectory of the growth index / has also been illustrated for the NGCG model with interaction. 
We find the value of / at a = 0.87 (i.e., z = 0.15) given by us agrees with the observational value of / = 0.51 ± 0.1 or 
/ = 0.58 ± 0.11 at the effective redshift of the survey z = 0.15. By comparing the NGCG and the GCG models with 
the ACDM model, we have shown that at early times, the growth indices of these three models arc the same. Near 
a — 0.09, the evolutionary trajectory of ACDM model would be different from ones of the other two. Near a ~ 0.2, 
the evolutionary trajectory of GCG model would deviate from one of the NGCG model. 
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